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Abstract. We discuss a method of constructing solution of the initial value problem for diffusion- 
type equations in terms of solutions of certain Riccati and Ermakov-type systems. A nonautonomous 



, Burgers-type equation is also considered. 

1. Introduction 

r-| ■ A goal of this note, complementary to our recent paper |37], is to elaborate on the Cauchy initial 
P^! value problem for a class of nonautonomous and inhomogeneous diffusion-type equations on R. A 
^ ' corresponding nonautonomous Burgers-type equation is also analyzed as a by-product. Here, we use 
' explicit transformations to the standard forms and emphasize natural relations with certain Riccati 
and Ermakov-type systems, which seem are missing in the available literature. Similar methods are 
' — ^, applied to the corresponding Schrodinger equation (see, for example, [B], [T], [S], [H], [21], 
CN ! ES], |3S], |SB], ISH] and references therein). A group theoretical approach to a similar class of 
^ partial differential equations is discussed in Refs. [12], [2H] and [51] . 

^ \ For an introduction to fundamental solutions for parabolic equations, see chapter one of the book 
^ by Friedman fH]. Among numerous applications, we only elaborate here on an important role of 
fundamental solutions in probability theory [lOj, f2T\. Consider an Ito diffusion X = {Xf '■ t > 0} 
O ■ which satisfies the stochastic differential equation 

dXt = b{Xt,t) dt + a{Xt,t) dWt, Xo = x, (1.1) 

in which W = {Wt : t > 0} is a standard Wiener process. The existence and uniqueness of solutions 
^ of (11. ip depends on the coefficients b and a. (See Ref. [21] for conditions of unique strong solution 
^ to (11.11) .) If the equation (II. ip has a unique solution, then the expectations 

u{x,t) = E, [0 (Ai)] =E[<f) (Xt) \Xo = x] (1.2) 

are solutions of the Cauchy problem 

Ut = ^(T^ {x,t)u^^ + b{x,t)u:,, u{x,0) = (t>{x) . (1.3) 

This PDF is known as Kolmogorov forward equation [TU], [2T]. Thus if p (x, y, t) is the appropriate 
fundamental solution of (II. 3p . then one can compute the given expectations according to 

E4(j){Xt)] = f p{x,y,t)(j){y) dy. (1.4) 
Jn 
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In this context, the fundamental solution is known as the probability transition density for the 
process and 

/ P{x,y,t) dy = 1. (1.5) 
Jn 

See also Refs. [1] and [20] for applications to stochastic differential equations related to Fokker- 
Planck and Burgers equations. 

2. Transformation to the Standard Form 

We present the following result. 
Lemma 1. The nonautonomous and inhomogeneous diffusion-type equation 

^ = « (^) ^ - (^) - c W ^) ^ + + f{t)x-h (t) x^) u, (2.1) 

where a, b, c, d, f, g are suitable functions of time t only, can be reduced to the standard autonomous 
form 

dr ~ de ^ ' 

with the help of the following substitution: 

u (x, t) = _^^e'''^tW+m-+<t)^ ^) ^ (2.3) 
V/^ if) 

i = P{t)x + e{t), T = ^{t). 
Here, fi, a, /3, 7, S, e, k are functions of t that satisfy 



^ + 2aa + d = (2.4) 



2/i 

and 

— + b-2ca- ^ao? = 0, (2.5) 
dt 

dB 

- {c + 4aa) (3 = 0, (2.6) 

^ - a/32 = (2.7) 

— -{c + Aaa)S = f -2ag, (2.8) 

(JiV 

de 

— + {g- 2a6) /3 = 0, (2.9) 

^ + gS-a5' = 0. (2.10) 

Equation (12. 5p is called the Riccati nonlinear differential equation [32], [ID], |12] and we shall 
refer to the system fl2.5l) - fl2.10l) as a Riccati-type system. 

The substitution (12.41) reduces the nonlinear Riccati equation (12. 5p to the second order linear 
equation 

/x" - r (t) /i' - 4a (t) /X = 0, (2.11) 
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where 



r (t) = — + 2c - Ad, a it) = ah + cd - + - (— - % 
a 2 \ a d 



(2.12) 



which shall be referred to as a characteristic equation [3( 

It is also known [37] that the diffusion-type equation (12. ip has a particular solution of the form 

(2.13) 



u 



^a(t)x'^+P{t)xy+-y(t)y'^+5(t)x+e{t)y+K{t) 



provided that the time dependent functions /i, a, 7, S, e, k satisfy the Riccati-type system (12. 4p - 

A group theoretical approach to a similar class of partial differential equations is discussed in 
Refs. [15], [28] and H 



3. Fundamental Solution 

By the superposition principle one can solve (formally) the Cauchy initial value problem for 
the diffusion-type equation (12. ip subject to initial data u (x, 0) = ip (x) on the entire real line 
— 00 < X < 00 in an integral form 



u (x, t) = Kq (x, y, t) (x) dy 



with the fundamental solution (heat kernel) 

1 



Kq (x, y, t) 



ao it)x^+l3o it)xy+"/o {t)y'^+So (i)x+£o it)y+Ko (*) 



V27r/io (t) 

where a particular solution of the Riccati-type system (I2.5p -( 12TT0|) is given by 

1 /x'o(t) d{t) 



7o it) 



ao (t) 

hit) 
/io it) ' 



4a it) jiQ it) 2a (t) ' 
h (t) = exp (^j^ (c (s) - 2d (s)) ds 



d (0) a (t) (t) _ p a{s)a{s)h{s) 



2a (0) /io(t)yt^'o(t) 7o (/i'o(s)) 
rf(0) 1 /ii(t) 



5o (t) 



2a (0) 2/Xi (0) /io (t) 

/ (^) + TTTV^? (^) /"o (^) + £Vt/^o (^) 



2a(t)/i(t) ^ 
^0 it) = ^77^ — ^0 it) - 



+ 2 



f^'o (t) 
^ a{s)h (s) 
t^'o is) 



a[s) - - y — - - 2a (s 

* a (s) 0" (s) h (s 
(^0 is)) 



ds 

2 (/io(s)(5o(s)) ds 



fis) + -T^gis) 

a [s) 



ds. 



(3.1) 
(3.2) 

(3.3) 
(3.4) 

(3.5) 
(3.6) 

(3.7) 
(3.8) 
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ft 



Kq {t) 



a jt) ^0 jt) 



^0 it) - 4 



a[s]a[s) / w / nn2 , 
(/^o {'V oo [s)) ds 







(3.9) 



+ 2 



as 



(/io(s) 5o(s)) 



as 



ds 



'0 /^ol^J 

with 5 (0) = g (0) / (2a (0)) , e (0) = —5 (0) , k (0) = 0. Here, /Iq and fii are the so-called standard 
solutions of the characteristic equation fl2.1ip subject to the following initial data 



/Xo(0) = 0, /i[,(0) = 2a(0)^0 /i^ (0) ^ 0, fi[{0) = 0. 



(3.10) 



Solution f i3.3p - fi3.9p shall be referred to as a fundamental solution of the Riccati-type system fl2.5p - 
f l2.10p : see fl3.27p - fl3.3ip and f l3.32p for the corresponding asymptotics. 



Lemma 2. The Riccati-type system ( [i^.^| )-( f^!T^) /ias the following (general) solution: 

/x(t) = -2/i(0)/io (t) (a (0) + 7o W) , 

/3oW 



a (t) = ao (t) 



4(a(0)+7o W)^ 
(5 (0) /3o (t) 



7(t)=7(0)- 
6{t) = 5^ (t)- 

£(t)=£(0)- 



2 (a (0)+7o it)y 
f (0) 



4 (a (0) + 7o (t)) 

/3o (t) (^(0)+eo (t)) 

2(a(0)+7o(t)) ^ 
/3(0) ((5(0) + 50 (t)) 

2(a(0)+7o(t)) ' 



K (t) = K (0) + (t) - 



(5(0)+5o(t))^ 



(3.11) 
(3.12) 

(3.13) 

(3.14) 

(3.15) 
(3.16) 
(3.17) 



4 (a (0)+7o {t)) 

in terms of the fundamental solution ^3. 3|) -( rOI) subject to arbitrary initial data ;U (0) , a (0) , /3 (0) 
7(0), 5(0), 6(0), /t(0). 

Proof. Use f l2.13p -f l3l2|) . uniqueness of the solution and the elementary integral: 



dy 



a 



a > 0. 



(3.18) 

Computational details are left to the reader. □ 
Remark 1. It is worth noting that our transformation l{2.3\) . combined with the standard heat kernel 



exp 



4(r-ro) 



(3.19) 



^Ati (r - To) 

for the diffusion equation h2.2) and I{3.11\) - f3.1'l\ ), allows one to derive the fundamental solution 
Ii3. 2\) of the diffusion-type equation Ii2.1\) from a new perspective. 
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Lemma 3. Solution ^3.11\) - ^.17\ ) implies: 



/So 



4(7-7(0)) 
(3(0) (3 



2(7-7(0))' 



7o = (0) - 



(0) 



and 



So = -5(0) + 



4(7-7(0)) 

/3(g-g(0)) 
2(7-7(0))' 

/3(0) (e-5(0)) 



which gives the following asymptotics 

1 

7o W 



2(7-7(0)) 
4(7-7(0)) 

c(0) a'(0) 



+ 



4a(0)t 4a (0) 8a2 (0) 
1 a' (0) 



+ 0(t) 



2a(0)t 4a2(0) 



0{t) 



1 c (0) a' (0) 



4a (0) t 4a (0) Sa^ (0) 



" ^ ^ 2a (0) 
/to (t) = (t) 



+ O (t) , So (t) 



2a (0) 



as t 0+. 



0(t) 



3.20) 
3.21) 
3.22) 
3.23) 

3.24) 
3.25) 
3.26) 

3.27) 

3.28) 

3.29) 

3.30) 
3.31) 



(The proof is left to the reader.) 

These formulas allows to establish a required asymptotic of the fundamental solution (13. 2p : 



Kq (x, y, t) 



X exp 



47ra (0) t 
a' (0) 



exp 



- yf 
Aa (0) t 



(3.32) 



_8a2 (0) 



^ y> 4a (0) 



exp 



_2a (0) 



{x-y) 



(Here, / ~ (7 as t — t- 0+, if limt_>o+ {f Iq) = 1- The proof is left to the reader. 



By a direct substitution one can verify that the right hand sides of f l3.lip - fl3.17p satisfy the 
Riccati-type system fl2.4p - fl2.10p and that the asymptotics f l3.27p -f l3.3ip result in the continuity 
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with respect to initial data: 

lim u(t) = u (0) , lim a (t) = a (0) , etc. (3.33) 

The transformation property fl3.11l) - fl3.17l) allows one to find solution of the initial value problem in 
terms of the fundamental solution fl3.3p - fl3.9l) and may be referred to as a nonlinear superposition 
principle for the Riccati-type system. 



4. ElGENFUNCTION EXPANSION AND ERMAKOV-TYPE SYSTEM 

With the help of transformation (12. 3p one can reduce the diffusion equation (12. ip to another 
convenient form 

which allows to find solution of the Cauchy initial value problem in terms of an eigenfunction 
expansion similar to the case of the corresponding Schrodinger in Refs. [25 ^ind [38] . This method 
requires an extension the Riccati-type system (I2.5p - (l2.10p to a more general Ermakov-type system 
[24j , which is integrable in quadratures once again in terms of solutions of the characteristic equation 
(12. lip . Further details are left to the reader. 

5. NONAUTONOMOUS BURGERS EQUATION 

The nonlinear equation 

dv f dv d'^v\ ( dv \ dv 

= 2(2b(t)x-f(t)), 

when a = 1 and b = c = f = g = 0,is known as Burgers' equation |2], [3], [5], [T7], [19], [35], [H] 
and we shall refer to (15. ip as a nonautonomous Burgers-type equation. 

Lemma 4. The following identity holds 

vt + a {vvr, - Vr^r,) + (fi- " cx) v^- cv + 2{f - 2bx) 

= -2 h - , (5.2) 



-2 — (The Cole-Hopf transformation) (5.3) 
u 



if 
and 

Qu = auxx — {g — cx) Ux + [d + fx — bx'^) u (5.4) 
(a, b, c, d, f, g are functions oft only). 

(This can be verified by a direct substitution.) 

The substitution (15. 3p turns the nonlinear Burgers-type equation (15. ip into the diffusion-type 
equation (12. ip . Then solution of the corresponding Cauchy initial value problem can be represented 

as 

(5.5) 



d 

V (x, t) = —2— In 
dx 



j Ko{x,y,t)exp (^-^ v{z,0) dz^ dy 
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where the heat kernel is given by fl3.2p . for suitable initial data v (z, 0) on M. 
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6. Traveling Wave Solutions of Burgers-type Equation 

Looking for solutions of our equation ( 15. ip in the form 

v = l3{t)F{l3{t)x + -i{t)) = l3F{z), z = f3x + j (6.1) 
(/3 and 7 are functions of t only), one gets 

F" = (co + ci) F' + FF' + 2c2Z + C3 (6.2) 

provided that 

(3' = cf3, i = (6.3) 

g = Cia/3, b = -^C2a(5^, (6.4) 

/ = ^0/3=^ (2c27 + C3) (6.5) 
(co, ci, C2, C3 are constants). From (16. 2p : 

F' = (co + ci) F + + C2Z^ + C3Z + C4, (6.6) 
where C4 is a constant of integration. The substitution 

F = -2— (6.7) 

transforms the Riccati equation (16. 6p into a special case of generalized equation of hypergeometric 
type: 

- (co + ci) A*' + 2 ('^2^^ + C32 + C4) yu = 0, (6.8) 



which can be solved in general by methods of Ref. [30]. Elementary solutions are discussed, for 
example, in [22] and |23] . 

7. Some Examples 

Now we consider from a united viewpoint several elementary diffusion and Burgers equations that 
are important in applications. 

Example 1 For the standard diffusion equation on M : 



dt dx"^ ' 

the heat kernel is given by 

\2 



K (x, t) = —= exp 



a = constant > (7.1) 
(x - vY 



Aat 



t > 0. (7.2) 



(See [1], [22] and references therein for a detailed investigation of the classical one-dimensional heat 
equation.) 
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Example 2 In mathematical description of the nerve cell a dendritic branch is typically mod- 
eled by using cylindrical cable equation [TS] : 

(7.3) 



du , 2 

r— = A 7— + u, T = constant > 0. 
at ox^ 



The fundamental solution on M is given by 

Kq (x, y, t) = ^— ^ exp 



t{x- vY 



t > 0. 



(See also (TB] and references therein.) 



(7.4) 



Example 3 The fundamental solution of the Fokker- Planck equation 



du d'^u ^ du ^ 
dt dx'^ dx 



on M is given by [37] : 

Kq (x, t) 

Here, 



v/27r(l -e-2*) 



exp 



lim Kq {x,y,t) = , 

t-i-co ^27C 



(x — e *?/)^ 
2(1 -e-2*) 



, t>0. 



y = constant. 



(7.5) 



(7.6) 



(7.7) 



Example 4 Equation 



du d'^u , , ^du 

~ = + yd - kx) a,k>0, g>0 



(7.8) 



dt dx"^ ' ^ dx' 

corresponds to the heat equation with linear drift when g = |28]. In stochastic differential 
equations this equation corresponds the Kolmogorov forward equation for the regular Ornstein- 
Uhlenbech process [TU]. The fundamental solution is given by 

Vke'^'/^ 



Kq (x, y, t) 



X exp 



A/47rasinh {kt) 
{k, {xe 



-kt/2 _ „,pkt/2 



+2^sinh(H/2))^ 



Aak sinh [kt) 



(See Refs. [10] and ^\ for more details.) 



t > 0. 



(7.9) 



Example 5 The viscous Burgers equation [2|, [3], [19], [23], [41] : 

a = constant > 



dv dv d'^v 
dt dx dx"^ ' 



can be linearized by the Cole-Hopf substitution [5] , [T7] : 

2a du 



(7.10) 
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which turns it into the diffusion equation (17. ip . Solution of the initial value problem has the form: 



V (x, t) 



a d 



In 



for t > and suitable initial data on 



exp 



Aat 



2a Jo 



V {z, 0) dz I dy 



(7.12) 



Example 6 Equation (I7.10p possesses a solution of the form: 

V = F {x + Vt) , V = constant 
(we follow the original Bateman paper [2J with slightly different notations), if 

VF' + FF' = aF", 

or 

(F + Vf ±A^ = 2aF', 
where A is a positive constant. The solution is thus either 

'A(x + Vt-c 



(7.13) 

(7.14) 
(7.15) 

(7.16) 

A + . + V - ^ ' ^'-''^ 

according as the + or — sign is taken. In the first case there is no definite value of v when a tends 
to zero, while in the second case the limiting value of v is either A — V or A + V according as x + Vt 
is less or greater than c. The limiting form of the solution is thus discontinuous [2] . 



or 



V + V = Atan 
A-v-V 



exp 



2a 

-(x + Vt 
a 



Further examples can be found in Refs. [TO], [23], [26], [28] and 



8. Conclusion 



In this note, we have discussed connections of certain nonautonomous and inhomogeneous diffusion- 
type equation and Burgers equation with solutions of the Riccati and Ermakov-type systems that 
seem are missing in the available literature. Traveling wave solutions of the Burgers-type equations 
are also discussed. 

Acknowledgments. We thank Professor Carlos Castillo- Chavez and Professor Carl Gardner for 
support, valuable discussions and encouragement. 
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